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ABSTRACT 


An  analysis  is  undertaken,  both  mathematical  and  numerical,  of  the 
quenching  of  two-dimensional  premixed  flames  under  a variety  of  circumstances. 
The  discussion  is  divided  into  two  parts,  the  first  of  which  deals  with 
quenching  due  to  proximity  to  a surface  through  which  there  are  heat 
losses,  the  second,  quenching  due  to  a shear  flow  of  the  kind  experienced 
by  a flame  attached  to  a wire.  The  discussion  includes  a critical  appraisal 
of  flame  stretch  and  the  role  sometimes  daisied  for  it  in  intuitive  explana- 
tions of  flame  quenching. 
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SIGNIFICANCE  AND  EXPLANATION 
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Premixed  flames  are  flames  that  occur  when  the  reactants  are  thoroughly 
mixed  before  combustion,  the  inner  cone  of  a bunsen  burner  flame  being  a good 
example.  They  are  distinct  from  diffusion  flames  (e.g.  a candle)  for  which 
the  fuel  and  oxygen  are  essentially  separated,  mixing  together  by  diffusion 
only  at  the  flame  itself.  It  is  possible  to  put  a flame  out  (quench  it) 
by  cooling  it  and  this  can  be  done  in  several  ways.  Inserting  a piece  of 
metal  (which  is  a good  conductor  of  heat)  into  the  flame  will  remove  heat 
from  it  for  example.  Wire  gauze  is  very  effective  in  this  respect  and  was 
used  for  this  very  purpose  in  the  miners  safety  lamp  developed  in  the  19th 
Century.  Blowing  on  a flame  in  the  right  way  can  also  cool  it. 

The  purpose  of  the  present  paper  is  to  examine  a mathematical  model  of 
a premixed  flame  and  investigate  the  precise  manner  in  which  the  two  cooling 
mechanisms  can  cause  quenching. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
suRinary  lies  with  MRC,  and  not  with  the  author  of  this  report. 
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Introduction 

The  equations  solved  in  this  paper  have  the  general  form 
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with  upstream  conditions 
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(1.1b) 


x-y-oo  T-^T,  Y-*l.  (1.2) 
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Additional  boundary  conditions  will  be  introduced  in  due  course.  Here,  T 
is  the  temperature,  Y the  mass  fraction  of  reactant,  L the  Lewis  Number, 

0 the  activation  energy,  and  Uf (y)  the  gas  velocity  which  is  assumed  to  be 
directed  in  the  positive  x direction  and  depends  only  on  y.  The  scaling 
is  such  that  when  L « 1 and  f(y)  =1,  a one-dimensional  flame  perpendicular 
to  the  X axis  has  a speed  (U)  of  unity  in  the  limit  0 -*•  «>. 

These  equations,  incorporating  as  they  do  the  fundamental  elements  of 
diffusion,  convection  and  reaction,  have  been  widely  used  in  the  study  of 
premixed  flames.  So  much  so,  that  they  have  acquired  a scientific  life  of 
their  own,  independent  in  some  respects  of  real  flames.  Thus  although  the 
problems  to  be  considered  in  the  present  paper  can  and  will  be  provided  with 
some  physical  motivation,  it  should  be  emphasized  that  our  primary  concern  is 
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with  an  understanding  of  the  nature  of  the  solutions  to  these  equations,  and 
not  with  the  explanation  of  experimentally  observed  phenomena.  In  this 
connection  we  eschew  the  all  too  common  practice  in  combustion  science  of 
viewing  the  relation  between  mathematics  and  experiment  as  a one  way  street  in 
which  the  only  role  of  mathematics  is  to  explain  known  experimental  facts. 

For  one  thing  mathematical  models  such  as  (1.1)  can  be  used  to  explore  the 
validity  of  intuitive  physical  arguments.  If  a mathematical  model  contains 
precisely  the  ingredients  upon  which  a physical  argument  is  based,  and  yet  has 
solutions  that  contradict  the  intuitive  reasoning,  then  the  latter  must  be 
abandoned.  One  of  the  concerns  of  this  paper  is  to  apply  such  a test  to  the 
concept  of  flame  stretch. 

One  of  the  ways  in  which  equations  (1.1)  are  deficient  as  models  is  that 
they  contain  no  fluid  mechanics.  It  is  always  conceivable  that  the  solutions 
of  such  a system  have  qualitative  features  quite  distinct  from  those  that 
would  be  found  if  more  accurate  equations  were  adopted,  but  this  possibility 
seems  unlikely  in  the  present  case.  The  techniques  of  the  present  paper  can 
be  extended  to  equations  that  correctly  incorporate  the  fluid  mechanics  so  that 
it  will  be  possible  to  examine  this  question  in  future  studies. 

Activation  Energy  Asymptotics  - A Dichotomy.  The  only  mathematical  tool  devised 
so  far  for  the  rational  solution  of  ecjuations  such  as  (1.1)  is  activation 
energy  asymptotics,  i.e.  the  construction  of  solutions  in  the  limit  0 
The  pioneering  work  in  premixed  flames  is  the  one-dimensional  steady  analysis 
of  Bush  and  Fendell  (1970) , and  more  recently  Sivashinsky  (1974a,  1974b,  1975) 
and  Buckmaster  (1977)  have  obtained  results  for  both  three-dimensional  and 
unsteady  flames.  The  structure  of  these  flames  is  characterized  by  a 
flame-sheet  of  thickness  to  which  all  the  reaction  is  confined,  and 
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moreover,  the  analysis  demands  that  the  temperature  in  the  region  behind  the 
flame  sheet  (where  Y = 0)  differs  by  only  an  amount  from  the  adiabatic 

value  (l+T^) • This  imposes  strong  restrictions  on  the  class  of  problems 
that  can  be  solved  in  this  fashion  and  it  is  useful  to  examine  these  limitations 
Defining  the  location  of  the  fleune  sheet  by 

X = x^(y)  (1.3) 

it  is  convenient  to  replace  x by  the  new  variable  s = (x-x^)  so  that  the 
flame-sheet  is  fixed  at  s = 0,  and  the  unburnt  mixture  occupies  the  region 
s < 0.  Adding  equations  (1.1)  eliminates  the  reaction  terms  and  we  may  write 

0+3  0+ 

Uf(y)  / (T+y)ds  = / (AT+^  AY)ds 


which  upon  evaluation  becomes 
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representing  a global  energy  balance. 

As  noted  above,  in  order  to  carry  out  a self-consistent  asymptotic  analysis 
each  of  the  terms  on  the  left  side  of  this  equation  can  have  a magnitude  of  at 
most  o(^)  » so  that  each  of  the  two  terms  on  the  right  must  be  similarly 
bounded.  Restricting  attention  to  situations  which  permit  large  flame  deforma- 
tions (i.e.  xj  « 0(1)),  there  are  two  ways  to  ensure  this.  The  first  is  to 
look  for  disturbances  that  are  described  on  an  0(6)  scale  by  writing 


x 


f 


(1.5) 


the  key  to  the  work  of  Buckmaster  (1977)  and  that  of  Sivaskinsky  referenced 


earlier.  The  second  is  to  unrite 


L = 1 . 0(i 


(1.6a) 


Y 


and  confine  attention  to  those  problems  for  which,  consistent  with  this, 


T + y 


1 + T 

oo 


(1.6b) 
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the  )cey  to  the  stability  analysis  of  Sivashinslcy  (1977)  . At  the  present  time 
a fair  amount  is  )<nown  about  solutions  generated  by  the  first  choice,  but  their 
role  in  combustion  theory  is  unclear  since  they  have  a propensity  for  local 
instability.  Much  less  is  )cnown  about  the  second  type  of  solution,  although 
the  stability  question  appears  more  favorable  (there  is  a band  of  Lewis  Numbers  for 
which  the  one-dimensional  flame  is  stable  for  this  class  of  disturbances) , and 
for  the  most  part  the  analysis  of  the  present  paper  is  confined  to  solutions 
of  this  type. 


* 

Though  not  clearly  stated  in  that  work. 


2 . Flame  Near  a Quenching  Surface 

The  effect  of  heat  losses  on  flames  is  of  great  practical  interest  and 

has  been  studied  for  many  years.  Sound  mathematical  results  are  restricted  to 

the  case  of  volumetric  heat  losses,  a problem  of  less  interest  than  one  in  which 

the  losses  occur  at  a boundary,  the  subject  of  this  section. 

The  problem  to  be  studied  can  be  thought  of  as  having  its  physical  origins 

in  the  quenching  of  a flame  at  a burner  rim.  The  flame  is  confined  to  the 

region  y > 0 (Fig.  1)  with  the  quenching  surface  located  at  y = 0.  This 

surface  is  assumed  to  be  chemically  inactive  so  that  the  diffusion  mass  flux  to 

it  is  zero,  but  g-j  heat  losses  are  permitted  at  a rate  proportional  to  the 

difference  between  the  surface  temperature  and  T , the  ambient  temperature 

00 

of  the  fresh  mixture.  For  large  values  of  (-x)  the  flcune  is  far  from  the 

surface  and  unaffected  by  it.  Then  for  a uniform  flow  (f=l)  the  flame  is  one- 

dimensional  and  inclined  at  an  angle  arcsin(U  to  the  horizontal.  As  x 

increases  the  flame  approaches  the  surface  and  is  effected  in  two  ways.  There 

is  a geometric  effect  present  even  in  the  absence  of  heat  losses,  which  occurs 

(again  considering  a uniform  flow)  because  the  one-dimensional  structure  is  not 

compatible  with  an  adiabatic  wall  condition  except  in  the  case  U = 1.  Secondly, 

heat  losses  will  tend  to  reduce  the  flame  temperature  (the  temperature 

immediately  behind  the  flame  sheet)  and  this  will  tend  to  slow  the  flame  down. 

Lewis  and  von  Elbe  (1961,  p.  213)  in  a discussion  of  burner  rims  argue  that 

quenching  occurs  when  the  heat  losses  reduce  the  flame  speed  to  zero,  so  that 

the  flame  of  Fig.  1 will  be  horizontal  at  the  quenching  point.  How  obvious 

this  is  depends  on  ones  intuitive  powers  of  course,  but  it  is  worth  noting  that 

it  is  in  sharp  contrast  with  results  for  one-dimensional  flames  with  volumetric 

heat  losses  (Buc)cmaster , 1976),  for  which  extinction  occurs  when  the  losses  are 

-1/2 

sufficient  to  reduce  the  flame  speed  to  a fraction  e of  the  adiabatic 
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value.  For  this  reason,  because  it  has  never  been  shown  that  there  are  solutions 


of  (1.1)  consistent  with  Lewis  and  von  Elbe's  physical  description,  and  because 
the  interplay  between  the  geometrical  and  heat  loss  mechanisms  has  never  been 
explored,  we  study  this  problem. 

Activation  Energy  Asymptotics.  The  analysis  of  equations  (1.1)  in  the  asymptotic 
limit  9+00  is  by  now  straightforward,  and  is  a generalization  of  Bush  and 
Fendell's  (1970)  analysis  of  the  one-dimensional  deflagration  wave.  Reaction 
is  confined  to  a thin  flame  sheet  (the  flame  slcetched  in  Fig.  1)  in  which  the 
gradients  of  T and  V change  rapidly,  and  both  within  the  flame  sheet  and 
elsewhere  solutions  are  sought  as  series  in  inverse  powers  of  0,  i.e. 

C C + C,  -I-  ol  — 

001  (e2 

In  view  of  (1.6)  it  is  appropriate  to  write 


i=l-|x,  X=  0(1)  , 


T + Y = ^ ~ (1+T  ) + r + O 
00  0 1 


\0' 


(2.1a) 


(2.1b) 


so  that  adding  equations  (1.1)  together  yields  an  equation  for  valid 

everywhere  outside  of  the  flame  sheet. 


Uf(y)  ~ XAT„  . (2.2) 

dX  1 0 

Behind  the  flame  sheet  Y vanishes  identically  (there  is  no  other  mechanism  to 

terminate  the  chemical  reaction)  and  T.,  » 1 + T , the  adiabatic  flame 

0 ® 

temperature.  Ahead  of  the  flame  the  reaction  is  frozen,  i.e. 

8T 

X < X-  Uf(y)  . (2.3) 

f dX  0 

Tq  is  continuous  across  the  flame  sheet  but  its  gradient  jumps,  the  magnitude 
of  the  jiimp  being  determined  from  an  analysis  of  the  flame  sheet  structure, 
which  we  now  consider.  Thus  within  the  flame  sheet. 
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(2.4a) 


T 

Y 
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1 , , 
+ T + - T (n) 

oo  0 

w (n)  +0 


1 


(2.4b) 

(2.4c) 


(y  = h(x)  defines  the  location  of  the  flame  sheet  and  is  a preferable 
description  to  one  using  since,  as  we  shall  see,  the  latter  is  not 

necessarily  single  valued) , whence 


d^ 


dn 


^ = -2'^(1+T  )"'*[l+h'^]'^a)  exp 
dn" 


(1+T  )' 


(2.5) 


The  quantity  (t+u)  is  a linear  function  of  n so  that  matching  with  the 
solution  on  each  side  of  the  flame  sheet  implies  that  'f ^ is  continuojs  across 
the  sheet  and  so  may  be  said  to  have  the  value  there,  whereupon 


T + U)  = 


If 


(2.6) 


Integrating  the  equation  for  t obtained  by  eliminating  oo  between  (2.5) 
and  (2.6)  then  leads,  in  the  usual  way,  to  the  result 

_ ^ 

8T  r ' 2l  ~ 2 

y " h - 0 — ^ [l+h  exp 

Since  Y vanishes  to  all  orders  immediately  behind  the  flame  sheet,  is 

the  perturbation  flame  temperature.  The  notion  that  an  increase  (decrease)  in 
flame  temperature  is  associated  with  an  increase  (decrease)  in  flame  speed  goes 
bacl<  to  Mallard  and  Le  Chatelier  (1883),  so  that  is  a quantity  of  great 

physical  interest. 

It  remains  to  observe  that  since  the  equation 

uf  (y)  ^ = (l-^A)A.^  + ^ AAT 

dX  D U 

is  valid  everywhere,  there  is  a jump  in  V<P ^ across  the  flame  sheet  whose 


If 


2(1+T  ) 


(2.7) 
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magnitude  is  defined  by 


6(V^^)  = -A6(Vt^)  . 


(2.8) 


Equations  (2.2)  and  (2.3)  must  be  solved  subject  to  the  conditions  (2.7), 
(2.8),  the  continuity  of  and  across  the  flame  sheet,  and  appropriate 

boundary  conditions  at  the  quenching  surface.  The  latter  are  ta)<en  to  be 

3Y  „ 9T  a 

y = 0 ~ ~ (T-T_)  (2.9a,b) 


9T 

g 

9y 


9i^ 


= 0, 


1 


= a(T^-T  ) , 

9y  0 “ 


(2.9c,d) 


where  a is  a prescribed  0(1)  constant.  Such  small  heat  losses  are  imposed 
on  the  problem  by  the  requirement  (1.6b),  but  small  though  they  are,  through 
their  direct  effect  on  they  can  substantially  influence  the  flame.  In 

reality  the  heat  losses  through  a surface  will  depend  on  the  temperature  field 
within  the  solid  phase,  but  by  adopting  (2.9)  this  complication  is  eliminated 
without  compromising  the  essential  physics. 

The  Slow  Flame  Approximation.  Activation  energy  asymptotics  has  reduced  the 
mathematical  problem  to  an  elliptic  one  involving  a free  boundary  (the  flame 
sheet) . Other  than  the  obvious  mathematical  difficulties  of  such  a problem 
there  is  a serious  conceptual  difficulty  which  arises  from  the  fact  that  the 
solution  at  any  point  depends  on  the  solution  everywhere.  It  is  therefore 
necessary  to  construct  the  solution  downstream  of  any  quenching  point  that  might 
exist,  where  the  flame  sheet  is  terminated,  and  yet  what  form  this  solution 
must  ta)ce  is  unknown  at  the  present  time.  A resolution  of  this  difficulty  is 
to  assume  that  the  gas  speed  is  much  greater  than  the  adiabatic  flame  speed,  i.e. 

U >>  1,  f = 0(1),  for  then  x derivatives  are  much  smaller  than  y derivatives, 
the  governing  equations  reduce  to  parabolic  type,  and  these  can  be  integrated 
by  marching  downstream  towards  the  quenching  point.  The  situation  is  in  some 
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ways  analogous  to  that  of  determining  the  location  of  the  separation  point  in 
a body  immersed  in  a viscous  flow.  Use  of  the  Navier-Stokes  equations,  which 
are  elliptic,  presents  enormous  difficulties  for  this  problem  since  the  nature 
of  the  flow  downstream  of  the  separation  point  is  poorly  understood.  If  the 
Reynolds'  Number  is  very  large,  however,  the  flow  can  be  described  by 
Prandtl's  boundary  layer  equations,  which  are  parabolic,  and  these  can  be 
integrated  towards  the  separation  point  (Schlichting , 1955  p.  131) . 

Formally  the  simplification  is  achieved  by  introducing  a new  independent 
variable 


X = -X 


(2.10) 


and  then,  in  the  limit  U the  equations  reduce  to 

3T 


X < X., 


|x  - X^l  > 0 


f (y) 


0 0 


3X 


3y 


2 ' 


3^1  32v>^  3\ 


3y 


3y 


Moreover  the  flame  sheet  condition  (2.7)  simplifies  to 

3T„ 


y ->•  h - 0 


3y 


exp 


If 


so  that  (2.8)  can  be  written  in  the  form 


\ 

f ' 

P 

= -X6 

0 

= A exp 

1 

i 3y  , 

2(1+T  ) J 


"if 


2(1+T  ) 


(2.11a) 


(2.11b) 


(2.12) 


(2.13) 
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Solution  for  a Uniform  Flow,  X = 0.  This  simplified  problem  can  be  solved 
numerically  for  arbitrary  choices  of  A and  f (y) , but  it  is  advantageous 
to  restrict  attention  to  a Lewis  Number  of  one  (X  = 0)  and  velocity  distribu- 
tions f for  which  can  be  explicitly  expressed  in  quadratures.  Lewis 

Number  effects  will  be  explored  in  a different  but  related  problem  in  §4. 

Consider  a uniform  flow  (f  = 1)  . Far  from  the  wall,  where  is  small, 

the  solution  is  characterized  by  a flat  flame  inclined  at  45“  to  the  horizontal. 
The  structure  in  this  far-field  limit  is 

x+y-h„ 


X < X, 


T„  = T + e 
0 “> 


= 0 


(2.14a) 


T„  = 1 + T , 
0 “> 


Xf  = -y  + ho  , 


= 0 


(2.14b) 

(2.14c) 


and  this  provides  the  initial  condition  for  the  general  problem.  The  constant 
ho  is  chosen  so  that  at  x = where  the  flame  is  a distance  ho  from  the 


wall,  the  influence  of  the  latter  is  small  enough  to  validate  (2.14) 
numerical  integration  can  then  proceed  downstream  from  x = 0> 


The 


The  value  of  at  the  flame  sheet  is 


"fi.  (x) 

^f 


1/2 


X (T  - T (s,0)  ] 

f , ” 0 

j ds  exp 


(X-s) 


1/2 


r_i  h^x.)] 

[ 4 (x-s)J 


(2.15) 


where  h is  the  distance  between  the  flame  sheet  and  the  wall  (h(0)  = ' 

and  in  this  way  the  problem  is  reduced  to  one  for  To  and  h alone.  It  is 
clear  that  the  wall  temperature  will  exceed  T^  so  that  is  non-positive 

everywhere  on  the  fl^une  sheet. 

Useful  insight  into  the  nature  of  the  solution  can  be  obtained  by  using 
an  approximate  technique  that  was  developed  to  a very  sophisticated  level  in  the 
study  of  viscous  boundary  layers  prior  to  the  advent  of  electronic  calculating 
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machines  (Schlichting,  1955,  p.  201) . Integrating  the  equation  for 
between  the  wall  and  the  flame  sheet  yields  an  exact  energy  balance,  namely 


_d_ 

dX 


h 


i 

0 


dy 


exp 


If 


2(1+T  ) 


+ h'  (X)  (1+T_^)  . 


(2.16) 


The  assumption  is  now  made  that  a fairly  crude  approximation  to  T^,  provided 
it  satisfies  the  boundary  conditions  at  the  wall  and  the  flcune  sheet,  can 
generate  a good  approximation  to  the  left  side  of  this  equation.  Thus 
approximating  T^  by  a simple  polynomial. 


2 .. 

y h 

1 Gxp 

If 

2h  " 2| 

, . 2 

2 (1+T  ) 

L*  OO  V 

and  substituting  into  (2.16)  yields  the  approximate  equation 

2 , 1 ^If 

- hh  ~ -1  --  2 • 

^ ^ (1+T  )^ 


(2.18) 


This  is  not  accurate  when  h is  large  because  (2.17)  fails  to  reflect  the 

important  fact  that  T^  will  only  differ  from  T in  the  0(1)  neighborhood 

0 ® 

of  the  flame  sheet,  but  it  should  be  a reasonable  approximation  when  h is  0(1) . 

When  the  wall  is  insulated  (a=0)  vanishes  identically  and  (2.18)  has 

the  solution 

h^  » 3(XQ-X)r  Xq  = constant  (2.19) 

so  that  the  flame  intersects  the  wall  at  x “ Xq*  That  this  phenomenon  is  not 
just  a creature  of  the  approximation  is  suggested  by  the  fact  that  a local  solu- 
tion of  the  exact  problem  can  be  found  corresponding  to  intersection  of  the  wall 
and  the  fleune  sheet.  For  if  it  is  assumed  tliat  when  (Xq~X)  is  small  the 
sheet  is  described  approximately  by 

h CCXq-x)^'^^  (2.20) 


( 
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then  in  the  neighborhood  of  Xq  it  is  appropriate  to  write 


0 


1 + T + (X--X)^^^F(i;) 
00  0 


which  is  a solution  of  (2.11a)  provided 

F” 


I CF-  + I F = 0 , 


(2.21) 


(2.22) 


c = y(Xo-x) 


-1/2 


The  boundary  conditions  (2.9c),  (2.12)  require  that 

F' (0)  = 0,  F' (C)  = 1 


(2.23a,b) 


so  that  the  appropriate  solution  of  (2.22)  is 

C 


d^e 


AC  / 

B 


where  B is  a solution  of  the  transcendental  equation 


^=/  dC 
® 0 


(e  -1) 


B > 0 


and  A is  chosen  to  ensure  (2.23b).  As  with  any  local  downstrecun  solution  of 
a parekbolic  problan  that  does  not  malce  use  of  the  initial  conditions,  there  are 
undetermined  constants  such  as  C. 

An  exact  local  solution  of  this  type  can  also  be  found  when  a is  not  zero, 
assuming  that  'fi  is  locally  a constant. 

When  the  wall  is  not  insulated  (cf^) , ifi  does  not  vanish  and  ' is 
negative  since  heat  losses  to  the  wall  cause  a steady  decline  in  the  flame 
temperature.  The  magnitude  of  at  any  point  depends  on  a and  it  is 

apparent  that  if  a is  large  enough  the  right  hand  side  of  (2.18)  may  well 
loecome  positive  at  some  value  of  x>  corresponding  to  a local  increase  in  h. 

Such  an  increase,  assuming  it  does  occur,  corresponds  to  a negative  flame  speed 
in  the  sense  that  relative  to  the  gas  the  flame  propagates  towards  the  burnt  region. 
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These  various  conjectures  have  been  tested  by  numerical  integration. 


The  equation  for  T^  can  be  written  in  the  form 

V 0 0 0 

. h - hh'z  — — = — 

3X  3z  .2 

oz 


(2.24) 


h(x) 

which  must  be  solved  in  the  domain  x > Of  O^z^l  subject  to  the  initial 
condition  (2.14a)  and  boundary  conditions 

3T_ 


z = 0 


3z 


= 0 , 


T„  = 1 + T , 
0 “ 


!!o 

3z 


= h exp 


2(1+T  ) 


(2.25a) 


(2.25b,c) 


Bac)cward  differencing  in  x yields  an  equation  for  T^^j^(z)  = 

which  contains  h,  , r h(x.  ,)  as  well  as  )cnown  data  at  X..  An  iterative 
3+1  3+1  3 

procedure  was  adopted  in  which  guessed,  the  equation  for  + 

integrated  to  satisfy  the  boundary  conditions  (2.25a,b),  and  adjusted 

until  (2.25c)  is  satisfied.  This  proved  to  be  rapidly  convergent.  Integration 
in  the  z direction  was  underta)ten  by  central  differencing  followed  by  inver- 
sion of  the  tridiagonal  matrix  using  a Choles)ci  decomposition. 

Results  for  the  flame  shape  are  shown  in  Fig.  2 for  different  values  of  a 

with  T fixed  (=  0.2) . For  sufficiently  small  values  of  a the  flame  height 
00 

decreases  with  a monotonically  decreasing  slope  until  intersection  with  the 
wall  occurs  in  a manner  consistent  with  (3.10).  Increasing  the  magnitude  of  o 
eventually  leads  to  the  appearance  of  an  inflexion  point  and  the  curve  displays 
a flattened  portion  of  increasing  length.  Ultimately  at  some  critical  value  of 
a(*  6.9),  h has  a local  minimum  denoted  by  Q,  and  just  downstream  of  Q h 
is  an  increasing  function  of  x* 
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Those  portions  of  the  curve  for  which  h is  an  increasing  function  of  x 


r* 

are  physically  unusual,  since  they  correspond  to  a negative  flame  speed.  Their 
peculiarity  is  seen  most  clearly  perhaps  by  noting  that  the  present  problem  is 
identical  to  that  of  a one-dimensional  unsteady  flame  travelling  down  a tube  ’ 

whose  end  is  sealed  by  a quenching  surface  at  y = 0,  with  x playing  the 
role  of  time.  A negative  flame  speed  corresponds  to  reflection  of  the  flame 
from  the  end  wall.  This  author  knows  of  no  experimental  reports  of  negative 
flame  speeds  so  that  a plausible  hypothesis,  adopted  throughout  this  paper,  is 
that  the  flame  is  quenched  at  Q.  The  overall  picture  is  then  consistent 
with  the  physical  description  of  Lewis  and  von  Elbe  (1961,  p.  216).  Note  that 
without  this  hypothesis  the  present  results  provide  no  explanation  of  quench- 
ing, for  the  numerical  computations  can  be  continued  downstream  of  Q with  no 
difficulty  until  ultimately  the  flame  intersects  the  wall. 

# 

The  region  between  the  quenching  point  and  the  wall  is  known  as  the  dead 
space  and  at  the  critical  value  of  a this  has  a minimum  thickness  of 
approximately  0.9. 

Solution  for  a Shear  Flow,  i = 0.  In  examining  the  behavior  of  a flame  near 
a stationary  wall  it  is  more  realistic  to  consider  a shear  flow,  rather  than  a 
uniform  flow,  so  that  the  problem  is  now  reconsidered  with  the  choice 

f(y)  = y . (2.26)  f 

This  alters  the  calculation  in  two  ways.  First,  the  value  of  at  the  flame 

sheet  is  now  given  by 


1 


(2.28) 


V,  ~ , ,1/2 

h ~ (hQ-2x) 


with  a temperature  profile 


Tq  ~ T + exp[y  - (hQ-2x)^^^) 


(2.29) 


valid  when 


2 1/2 

X y - (h^-2x)  = 0(1)  . 


(2.30) 


For  as  x (y  fixed)  , -*■  T as  required;  when  y = h,  T„  = 1 + T 

0 0®  0 


3T. 


2 

and  -r-^  =1  as  required;  and  (2.29)  is  a solution  of  (h^-2x)^^^  — ^ ^ / 

3y  0 3X  . 2 

dy 


a valid  approximation  to  the  true  equation  for  in  the  region  defined 

by  (2.30). 

The  results  of  numerical  integrations  starting  with  the  initial  (x=0) 

profile  (2.29)  are  shown  in  Fig.  3.  Curves  are  drawn  for  different  values  of 

a with  once  again  fixed  at  0.2,  and  the  results  are  qualitatively  similar 

to  those  for  a uniform  flow.  The  local  structure  near  the  intersection  point 

of  flcune  and  wall  is  now  described  (when  = 0)  by 

^f 


h ~ C(Xq-x) 


1/3 


,1/3, 


where 


’’o  ^ (Xq'X)  ' F(C) 


F"  - j CF*  + J F - 0, 


(2.31a) 

(2.31b) 


(Xq-X) 


1/3 


F'  (0)  - 0,  P' (C)  - 1 . 
The  general  solution  for  P is 


F - AC  / dC 
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2.8 


FIG.  3.  EFFECT  OF  HEAT  LOSSES  FOR  A SHEAR  FLOW 


and  B must  be  chosen  so  that 


Z=/  dc 

B 0 


(e®  -1) 


B > 0 . 


Concluding  Remarks.  The  analysis  of  this  section  is  founded  on  the  assumption 
that  the  gas  speed  is  much  greater  than  the  flame  speed  so  that  y derivatives  are 
much  larger  than  x derivatives.  Clearly  this  condition  is  violated  in  the 
vicinity  of  the  intersection  of  wall  and  flame,  so  that  it  is  not  clear  how 
the  nonquenched  solutions  should  be  interpreted.  Certainly  for  the  shear  flow 
one  would  expect  flashback  to  occur  if  the  heat  losses  are  insufficient  to 
cause  quenching,  and  then  there  can  be  no  steady  solution.  These  considerations 
do  not  cast  any  doubt  on  the  prediction  of  quenching  of  course,  since  the 
quenched  solutions  have  bounded  flame  slopes  and  the  slow  flame  approximation 
is  uniformly  valid  all  the  way  to  the  quenching  point. 

Finally,  it  can  be  immediately  recognized  that  for  unquenched  flames  in  a 
uniform  flow  the  flame  speed  increases  as  the  wall  is  approached  despite  the 
fact  that  the  perturbation  flame  temperature  is  negative.  This  provides 

a counterexample  to  any  speculation  that  even  for  two-dimensional  flames  a 
decrease  (increase)  in  flame  temperature  can  always  be  associated  with  a 
decrease  (increase)  in  flame  speed. 


I 

I 
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3.  Some  Remarks  on  Flame  Stretch 


In  §2  it  is  shown  that  if  a flame  comes  close  to  a surface  through  which 
there  are  large  enough  heat  losses,  it  will  be  quenched.  The  fact  that  the 
surface  is  a phase  boundary  plays  no  role  and  it  may  be  conjectured  that  under 
some  circumstances  the  heat  transfer  from  one  part  of  the  combustion  field 
to  another  can  cause  quenching  even  in  the  absence  of  external  heat  sinks. 

Indeed  it  has  been  argued  that  an  appropriate  measure  of  this  heat  flux  is 
the  local  flame  stretch,  and  that  if  this  is  greater  than  some  critical  value 
quenching  will  occur. 

Before  going  any  further  it  is  appropriate  to  define  what  is  meant  by 
flame  stretch  and  to  do  this  in  an  uncunbiguous  fashion  it  is  necessary  to 
define  the  concept  of  a flame  surface.  This  is  simply  a sheet  that  characterizes 
the  location  of  the  flciroe.  Thus  for  the  problem  of  §2  the  reaction  flame 
sheet  is  a flcune  surface  when  looked  at  on  a scale  that  is  large  compared  to 
the  preheat  zone  thickness  divided  by  0.  For  similar  problems  Lewis  and 
von  Elbe  consider  the  flame  surface  defined  by  the  locus  of  inflection  points 
of  the  temperature  distribution  in  the  preheat  zone.  For  problems  of  the  type 
characterized  by  (1.5)  which  this  author  has  called  'slowly  varying  flames' 
(Buckmaster , 1977),  hydrodynamic  disturbances  are  described  on  a scale  that 
is  0(0)  times  the  preheat  zone  thickness,  and  on  this  scale  the  flame  is 
simply  a hydrodynamic  discontinuity  of  the  kind  discussed  by  Markstein  (1964)  . 
This  discontinuity  is  a flame  surface.  Now  consider  a point  moving  in  a flame 
surface  with  a tangential  velocity  equal  to  the  tangential  component  of  the 
gas  velocity  immediately  ahead  of  the  surface.  A small  area  A comprised 
of  such  points  will  in  general  be  deformed  as  it  moves  over  the  surface,  its 
magnitude  will  change  and  the  flame  will  be  stretched.  A precise  measure  of 

this  stretch  is  (Williams,  1975). 

A dt 
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Karlovitz  (1953)  were  the  first  to  suggest  that  this  stretching  can, 

if  large  enough,  lead  to  extinction  of  the  flame,  and  the  idea  was  subsequently 

* 

pursued  with  great  enthusiasm  by  Lewis  and  von  Elbe  (1961)  . They  argue,  for 

example,  that  this  is  the  cause  of  blowoff  of  a flame  stabilized  behind  a 
straight  wire  (Fig.  4).  Such  a flame  experiences  strong  shearing  from  the 
boundary  layer  generated  by  the  wire,  and  associated  with  this  shearing  is  a 
positive  degree  of  stretch.  The  greater  the  speed  of  the  gas  past  the  wire  the 
greater  the  shear  and  therefore  the  stretch,  and  the  claim  is  made  that  through 
this  mechanism  there  is  a maximum  gas  speed  beyond  which  the  flame  will  not  stay 
attached . 

Now  there  are  several  a priori  objections  that  may  be  raised  to  the 
various  claims  that  have  been  made  for  stretch,  no  less  important  because  they 
are  clearer  a posteriori.  The  English  language  is  an  imprecise  tool  so  that 
the  manner  in  which  stretch  achieves  what  is  claimed  for  it  is  not  too  clearly 
described.  It  will  suffice  to  say  that  the  basic  idea  expressed  on  p.  227  of 
Lewis  and  von  Elbe  is  that  stretch  changes  the  rate  at  which  heat  is  transferred 
from  the  reaction  zone  to  the  unburned  gas,  and  this  effects  the  flcune  speed. 
This  can  not  be  a local  phenomenon  in  general,  since  the  thermal  history  of 
fluid  that  travels  from  far  upstream  to  the  flame  sheet  must  depend  on  the 
stretch  experienced  by  a finite  portion  of  the  flame,  if  it  depends  upon  the 
stretch  at  all.  Thus  a quenching  criterion  expressed  in  terms  of  a local 
stretch  factor  ^ ^ can  only  be  meaningful  if  the  local  stretch  is  repre- 
sentative of  the  stretch  experienced  by  the  fleune  globally.  But  there  is  no 


* 

That  they  attach  great  significance  to  the  concept  is  suggested  by  the  fact 
that  it  is  mentioned  in  the  Preface  to  the  Second  Edition  of  their  boo)(. 
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FIG.  4.  FLAME  STABILIZED  BEHIND  A STRAIGHT  WIRE 
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way  of  relating  the  global  stretch  experienced  by  the  flame  of  Fig.  4 at  one 


value  of  the  gas  speed  and  wire  diameter  to  that  experienced  under  different 

* circumstances  since  the  flows  will  not  be  similar  in  qeneral.  Thus  if  it  was 
true  that  blowoff  occurs  once  the  global  stretch  experienced  by  the  flame 
exceeds  some  critical  value,  then  a series  of  experiments  on  wires  of  different 
dicuneters  with  the  gas  composition  fixed  would  inevitably  yield  different 
values  of  some  locally  defined  Karlovitz  Number.  Only  if  the  preheat  zone 
thickness  is  much  smaller  than  the  radius  of  curvature  of  the  flame  will  the 
effect  be  purely  local. 

Unfortunately  there  is  little  reason  to  believe  that  the  heat  transference 
rate  is  only  associated  with  the  stretch,  even  in  some  global  fashion. 

Variations  in  the  preheat  zone  thickness  will  also  play  a role  and  these  are, 
in  some  sense,  independent  of  stretch.  One  only  has  to  think  of  a one- 
. dimensional  unsteady  flame  for  which  the  stretch  is  identically  zero  and  yet 

variations  in  thickness  can  occur  affecting  the  heat  flux. 

Also  the  argument  on  p.  227  of  Lewis  and  von  Elbe  that  positive  stretch 
must  reduce  the  flame  speed  is  based  only  on  considerations  of  transfer  of  heat 
and  active  species  and  fails  to  account  for  the  fact  that  the  diffusion  of 
fresh  mixture  is  also  affected  by  stretch  but  that  this  flux  is  in  the  opposite 
direction  to  the  heat  flux.  The  influence  of  positive  stretch  through  this 
mechanism  is  apparently  such  as  to  increase  the  flame  speed  so  that  the  overall 
effect  is  Lewis  Number  dependent.  At  a closed  axisymmetric  flame  tip  the 

• stretch  is  negative  and  the  flame  speed  enhanced  over  the  adiabatic  value. 

But  at  an  open  axisymmetric  tip  the  stretch  is  still  negative  but  here  the 

flame  speed  is  below  the  adiabatic  value,  something  not  admitted  by  the  claim  j 

on  p.  227,  Curiously  enough  on  p.  218  Lewis  and  von  Elbe  admit  the  role  of 
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reactant  diffusion  effects  in  explaining  open  tip  flames,  but  mention  it 
nowhere  else  in  their  extensive  discussion  of  stretch.  Thus  on  p.  242  where 
they  consider  blow-off  from  a wire  they  argue  that  the  stretch  is  positive, 
the  flame  speed  reduced  as  a consequence,  and  for  large  enough  stretch  the 
flame  will  be  quenched.  To  be  consistent  with  their  discussion  on  p.  218  they 
would  surely  have  to  admit  that  for  some  mixture  compositions  blow-off  is 
impossible,  but  this  important  test  of  the  validity  of  the  blow-off  argument  is 
nowhere  discussed.  Moreover,  on  p.  337  where  ignition  by  a point  source  is 
discussed,  the  admission  that  for  some  mixtures  stretch  causes  an  increase  in 
the  flame  speed  would  lead  to  the  curious  conclusion  that  for  such  mixtures  an 
arbitrarily  small  ignition  source  will  always  lead  to  ignition. 

These  remarks  are  illuminated  by  a general  equation  governing  slowly  vary- 
ing flames  that  is  derived  in  the  Appendix.  This  equation  ((A. 19)  with  t 
replaced  by  t')  is 


= 6 


dV 

dt’ 


dA 

dt’ 


dH 

dt 


(3.1) 


where  H is  the  flame  speed  divided  by  the  adiedaatic  flame  speed,  6 is  a 
Lewis  Number  dependent  constant  defined  by  (A. 3b),  t'  is  physical  time,  and 
V is  the  volume  of  a flame  element  defined  as  the  product  of  A and  an 
appropriate  flame  thickness  6 which  emerges  in  a natural  way  from  the  large 
activation  energy  analysis.  This  relation  is  a local  one  simply  because  an 
essential  feature  of  slowly  varying  flames  is  that  the  thickness  of  the  preheat 
zone  is  much  smaller  than  its  radius  of  curvature. 

Equation  (3.1)  reveals  kxath  the  thickness  and  Lewis  Number  effects 
discussed  above.  Thus  it  is  not  possible  to  say,  without  cjualif ication,  that 
positive  stretch  will  always  slow  the  flame.  Nor  can  it  be  argued  that  when 
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S is  negative  (L>1)  positive  stretch  necessarily  slows  the  flame,  although 
that  is  certainly  true  when  the  dilatation  is  positive.  Thus  physical  arguments 
b.jsed  on  stretch  alone  omit  an  essential  part  of  the  physics  and  so  are 
unacceptable. 

It  is  worth  noting  in  passing  that  equation  (3.1)  suggests  the  possibility 

of  a quenching  mechanism  quite  different  from  that  described  in  §2.  Figure  5 

shows  how  changes  with  6 ^ it  is  apparent  that  the  latter  has  a 

V dt 

-1  -1/2 

minimum  value  of  (-e  ) when  H = e . Thus  no  flame  can  exist  for  a 


dilatation  that  satisfies  the  inequality 


-B  i ^ i , 

V dt  e 


(3.2) 


and  at  the  critical  value,  when  the  flame  is  just  quenched,  the  flame  speed  is 
not  zero.  Such  a quenching  mechanism,  assuming  it  exists,  is  analogous  to 
quenching  by  volumetric  heat  losses  (Buc)cmaster , 1976)  and  not  the  quenching 
mechanism  of  §2.  The  qualification  is  necessary  in  view  of  the  fact  that  the 
dilatation  is  not  something  that  can  be  directly  controlled  but  is  determined 
in  part  by  the  instantaneous  flame  configuration.  Thus  it  is  conceivable  that 
no  matter  what  shear  flows  or  other  disturbances  are  applied  the  flame  always 
adjusts  so  that  the  dilatation  is  smaller  than  the  critical  value.  No  flame 
solutions  have  yet  been  obtained  that  exhibit  quenching  of  this  type,  and  the 
question  is  presently  an  open  one.  True,  Sivashinslty  (1976)  has  proposed  that 
a slowly  varying  flame  located  in  a stagnation  point  flow  will  not  exist  under 
some  circumstances,  and  there  are  obvious  similarities  between  his  analysis  and 
the  quenching  mechanism  proposed  here,  but  his  discussion  is  flawed  by  unresolved 
nonuniqueness . 

So  far  our  discussion  of  stretch  has  been  restricted  to  broad  qualitative 
considerations.  The  quantitative  manner  in  which  stretch  has  been  used  to 
correlate  blowoff  data  exposes  another  flaw.  In  blowoff  situations,  stretch 
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VARIATION  OF  FLAME  SPEED  WITH  DIIATATION  RATE  FOR  SLOWI.Y  VARYING  FIAMES 


occurs  because  of  shear  and  it  is  argued  that  the  stretch  factor  — — 

A ri  t 

(which  is  difficult  to  determine  experimentally)  can  be  character ize--*  by  the 


Karlovitz  Number 


- il  dU’ 

“^a  ^ U'  dy' 


(3.3) 


where  6'  is  the  thickness  of  the  preheat  zone  for  the  one-dimensional 

adiabatic  flame,  U*  is  a characteristic  velocity  and  is  characteristic 

dy 

shear.  It  is  then  argued  (see  Reed  1967  for  example)  that  for  identical  burner 
configurations  blowoff  of  different  mixtures  will  occur  at  the  same  value  of 
K . What  is  curious  about  this  claim  is  that  the  only  role  the  Lewis  Number 
plays  is  an  indirect  one  through  its  effect  on  the  adiabatic  flame  speed  and 
therefore  6',  and  yet  one  would  surely  expect  the  solutions  of  equations  such 
as  those  discussed  in  §2  to  depend  directly  on  the  Lewis  Number  through  their 
dependence  on  A.  This  expectation  will  be  verified  in  §4  where  it  is  shown 
that  the  response  of  a flame  to  a shear  flow  can  depend  critically  on  the 
value  of  A . 

In  view  of  the  preceding  discussion  it  is  this  author's  firm  opinion  that 
use  of  a local  stretch  criterion  to  explain  or  correlate  quenching  is  at  best 
misleading,  and  at  worst  simply  wrong.  We  note  that  Melvin  and  Moss  (1973), 
on  the  basis  of  a careful  examination  of  experimental  data,  have  also  questioned 
the  concept  as  applied  to  the  understanding  of  blow-off  from  burner  ports. 

These  objections  leave  unanswered  the  interesting  question  of  wliether  or 
not  a shear  flow  such  as  that  of  Fig.  4 can  quench  a flame,  and  this  is  the 
subject  of  14. 


* 

However,  the  differences  between  K and  the  stretch  factor  are  revealed  in 
the  Appendix. 
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t^iuenching  by  a Shear  Flow 


In  this  section  we  return  to  the  question  raised  at  the  beginning  of  §3, 
n.amely  whether  or  not  distortion  of  a flame  can  generate  sufficient  heat 
transfer  from  one  part  of  the  combustion  field  to  another  to  cause  quenching. 

The  specific  problem  to  be  examined  is  partly  motivated  by  the  physical  situation 
sketched  in  Fig.  4 but  it  must  be  emphasized  that  the  results  of  this  section, 
despite  their  undoubted  physical  interest,  probably  have  no  relevance  to  the 
question  of  blowoff  from  a wire.  Undoubtedly  this  is  better  understood  in 
terms  of  the  destabilizing  influence  of  the  velocity  gradients  behind  the  wire, 
and  the  stabilizing  influence  of  heat  transfer  from  the  flcime  to  the  wire. 

The  model  problem  is  shown  in  Fig.  6.  A combustible  mixture  flows  from 
left  to  right  with  a uniform  velocity  in  the  region  y > 0 and  a linear  shear 
in  the  region  y < 0.  Far  upstream  the  shear  does  not  influence  the  flame 
which  is  then  one-dimensional  and  inclined  to  the  horizontal  at  an  angle  defined 
by  the  gas  speed  and  the  adiabatic  flame  speed.  As  x increases  the  flame 
approaches  the  shear  region,  becomes  curved,  and  the  combustion  field  is  two- 
dimensional.  The  fundamental  concern  is  whether  or  not  there  are  values  of 
shear  and  Lewis  Number  for  which  the  flame  will  reverse  and  move  away  from  the 
X axis  with  increasing  x,  a phenomenon  that  we  identify  with  quenching. 

The  basic  equations  are  identical  to  those  of  §2  after  the  large  activation 
energy,  slow  flame  approximations  are  adopted,  with  the  function  f(y)  which 
defines  the  variations  in  gas  speed  chosen  as  follows, 

f - 1,  y > 0 (4.1a) 

f » 1 - Ky,  y < 0 . (4.1b) 

Numerical  Integration.  The  numerical  problem  is  slightly  more  complicated  than 
that  of  §2  since  the  analysis  is  not  restricted  to  X ■ 0.  The  variable  y is 
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replaced  by  z H y - h(x)  in  order  to  fix  the  flame  sheet  at  z = 0,  and  the 
difference  equations  are  solved  using  a double  iterative  scheme.  In  this 
scheme  <p  is  guessed  and  the  problem  for  T defined  by  (2.11a),  (2.25b) 

If  0 

solved,  iterating  on  h until  (2.12)  is  satisfied.  Then  solving  for 
everywhere  ((2.11b),  at  the  flame  sheet)  the  condition  (2.13)  is 

checl^ed  and  if  not  satisfied  to  the  required  accuracy  a new  guess  for  ^ is 
adopted  and  the  entire  procedure  repeated.  This  wor)cs  well  enough  for  moderate 
values  of  A.  The  initial  profile  with  which  the  calculations  are  started  at 
X = 0 is 

z < 0 T„  = T + e^  (4.2a) 

= -Az  e^  (4 . 2b) 

z > 0 T„  = 1 + T (4.2c) 

0 ® 

= 0 . (4. 2d) 

Some  results  of  such  calculations  are  shown  in  Fig.  7 for  different  values 
of  A with  K fixed  at  the  value  5.  It  is  apparent  that  both  quenched  and 
unquenched  profiles  are  possible,  depending  on  the  Lewis  Number.  Thus  A = 10 
gives  rise  to  quenching,  but  A = 0,5  apparently  does  not.  For  some  values  of  K 
quenching  can  also  occur  for  negative  values  of  A,  as  we  shall  see,  but  no 
computer  solutions  have  been  obtained  for  such  cases.  These  solutions  are 
characterized  by  large  values  of  K and  (-A),  with  quenching  occurring  at  large 
values  of  x»  and  convergence  difficulties  experienced  in  an  attempt  to  find 
such  solutions  were  so  serious  that  the  attempt  was  abandoned. 

No  attempt  has  been  made  to  map  quenching  boundaries  in  the  A - K para- 
meter plane,  an  expensive  procedure  not  of  obvious  value.  Instead,  the  asymptotic 
limit  K -►  ■>  has  been  explored  to  lay  bare  some  of  the  complexities  of  the  response. 
Aaymptotic  Analysis  in  the  Limit  K ♦ ■».  Limit  equations  valid  as  K -►  “>  in 
y < 0 are 
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FIG-  7.  FLAME  IN  A SHEAR  FLOW 
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(4.3) 
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1 
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Thus  provided  the  flame  remains  in  the  region  y > 0 (i.e.  h > 0)  it  is 

appropriate  to  solve  the  equations  in  the  region  y > 0 (where  f = 1) 

subject  to  boundary  conditions 

= T , 'fi=0  at  y = 0 . 

0 “>  1 

This  system  has  a simple  x~independent  solution  that  can  be  expected  to  be 
appropriate  as  X namely 

X 


(4.4) 


h = h = exp 


y < h 


-X  exp 


2(1+T  ) 

"O 

-X 


'a,  - 


T„  = T + exp 

0 oo  ^ 


2(1+T  ) 
00 

-X 


2(1+T  )' 


(4.5a,b) 


(4.5c) 


(4.5d) 


y > h 


-X,  Tq  = 1 + T„  . 


(4.5e,f) 


Now  the  flame  is  first  affected  by  the  shear  flow  when  the  fringe  of  the  preheat 
zone  defined  by  (4.2)  first  intersects  the  x axis.  That  is  why  in  Fig.  7 h is 

less  than  2 before  significant  curvature  is  apparent.  On  the  other  hand  quite 

* 

moderate  values  of  positive  X can  yield  large  values  of  h , for  example 

* 

X = 6,  T^  = 0.2  gives  h * 8.  Thus  if  X is  large  enough  it  can  be  expected 
that  the  limiting  solution  valid  as  K -*■  «>  will  display  quenching.  Indeed 
numerical  results  show  that  quenching  does  occur  if  X > 0,  but  not  if  X £ 0 
(Fig.  8)  . 

The  limiting  solution  characterized  by  (4.4)  is  only  valid  when  X = 0(1), 
for  when  x “ 0(K)  there  is  significant  diffusion  in  the  region  y < 0. 

Limit  equations  valid  in  this  region  are  obtained  by  first  writing 
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FIG.  8.  FLAME  SHAPE  WHEN 


so  that  in  the  limit,  assuming  h is  positive  which  is  certainly  true  in  the 


neighborhood  of  o = 0,  we  have: 


y < 0 


-y 


3^0  3 ^0 


3o  . 2 ' 
3y 


d<fi 

1 1 . ^ v 

IT  ■ TT  * ‘ 


(4.7a,b) 


ay 


ay 


a^T„ 


0 < y < h 


0, 


= 0 , 


(4.7c,d) 


ay 


ay 


h < y 


= 1 + T , 
0 " 


= 0 


(4.7e,f) 


ay 


In  the  region  y > 0,  both  and  are  linear  functions  of  y characterized 


by  h and  ^ , 

^f 


h < y 


(4.8a) 


0 < y < h 


T„  = exp 
0 


2(1+T  ) , 


(y-h)  + 1 + T 


(4.8b) 


r 


= -X  exp 


2(1+T  )' 


(y-h)  + sfi 


(4.8c) 


On  the  other  hand,  in  y < 0 equations  (4.7a,b)  have  similarity  solutions 


which  satisfy  the  conditions  T -*•  T^,  ->•  0 as  y -*•  neunely 

1 0 oo  1 


T - T + o f(n)  , 
0 “> 


(4.9a) 


1/3  ' 


f - f(0) 


1 3 


dn 


9 


(1-e  ) 


(4.9b,c) 


“ 0 g(n)  , 


(4.9d) 


1 3 


fia\  ^9^  r ^ 

g - A q~  n / <3n  ne  + g(0)|  1 + n / 


1 3 


In 


9^ 


(l-e  ) 


(4.9e) 


It  follows  that 


f (0) 
f(0) 


q'  (0) 
g(0) 


= 9 


3,(1 


(4.10a) 


(4.10b) 


The  solutions  obtained  in  the  two  regions  must  match  in  the  sense  that  'P 

and  their  first  derivatives  with  respect  to  y are  all  continuous  at  y = 0. 

Thus  (4.8)  defines  f,  g and  their  first  derivatives  at  n = 0 in  terms  of  h 

and  whence  (4.10)  may  be  rewritten  in  the  form 

f 


Aa^  = 9 ^ r( jj  (1-hA) 

1 _ 5 1 

-AAa^  = -X9  ^rj-|j  (1-hA)  + 9 ^ r(||  (VJ^^  +hXA) 


exp 


2<1+T  ) , 

. * j 


(4.11a) 


(4.11b) 


(4.11c) 


These  equations  determine  how  both  h and  P vary  with  a and  it  is  easily 

verified  that  matching  with  (4.5)  is  achieved  as  a 0. 

Since  it  has  already  been  concluded  that  when  X is  positive  quenching 

occurs  on  the  x ° 0(1)  scale,  we  shall  only  discuss  the  solution  of  these 

equations  for  negative  values  of  X.  By  eliminating  h a single  equation  for 

p may  be  written  down 
^f 


Mil 


1 

r +x) 

r X 1 

3 

0 = (^,  +X)  exp 

^f 

1 2 

2 

L 2(1+TJ^J 

^f 

2{1+T  ) 

_ 00  ^ 

(4.12) 


whence  it  follows  that  >P  decreases  monotonically  from  its  value  (-X)  at 
a * 0,  approaching  (-«“)  as  a “ (although  it  must  be  remembered  that  the 
analysis  is  only  valid  for  the  interval  in  which  h is  positive) . Turning  to  h. 


-35- 


(4.11)  yields 


- - Ah'  . 


"i,  r 

2(l-t-T 


2(1+T  ) ^ ^ 

“ 3 ] 


(4.13) 


from  which  it  may  be  concluded  that  if  A satisfies  the  inequality 


A < -6(1+T  )' 


(4.14) 


h is  an  increasing  function  in  the  neighborhood  of  a = 0 so  that  quenching 

* 

occurs  with  a dead  space  equal  to  h . 

2 

When  A lies  in  the  interval  (-6(1+T  ) ,0]  h decreases  from  its  value 

00 

* 

h at  o = 0 and  vanishes  when  o = o . where 

i 


'^i  = itlfl]  “ 


(4. IS) 


(1+T  )' 


a result  that  can  be  chec)ced  numerically.  Thus  in  Fig.  9 values  of  x(h=0) 
for  different  values  of  K are  shown  for  A = 0,  -1  and  through  these  points 
are  drawn  straight  lines  of  slope  equal  to  values  of  o greater  than 

the  flcune  sheet  lies  in  the  region  y < 0 and  its  progress  must  be  followed 
numerically  since  the  simple  similarity  solutions  (4.9)  are  no  longer  appropriate. 
Typical  results  are  shown  in  Fig.  10  and  show  no  evidence  of  quenching. 

Concluding  Remar)ts.  The  results  of  numerical  calculations  such  as  those  shown 
in  Fig.  7,  and  the  asymptotic  analysis  when  K “ suggest  that  for  some  range 
of  values  of  K there  is  a band  of  Lewis  Numbers  for  which  quenching  does 
not  occur,  but  outside  this  band  the  flcune  is  quenched  by  the  shear  flow.  The 
evidence  suggests  that  the  band  always  contains  the  value  A = 0.  This  appears 
to  be  the  first  sound  theoretical  evidence  that  a shear  flow  can  quench  a flame. 

Quenching  when  A is  positive  is  characterized  by  a monotonic  decrease  in 
the  flame  temperature  between  x • 0 and  the  quenching  point.  However,  if 


CONFIRMATION  OF  THE  ASYMPTOTIC  RESULTS  IN  THE  LIMIT 


quenching  occurs  when  X is  negative  first  increases,  reaching  some 

positive  maximum  value,  and  then  decreases  becoming  negative  before  the 
quenching  point  is  reached.  There  appears  to  be  little  doubt  that  negative 
values  of  are  necessary  for  quenching  to  occur. 
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—■  Sununary 


In  this  paper  the  explicit  analysis  of  two  problems  and  the  reinterpreta- 
tion of  certain  results  for  slowly  varying  flcunes  has  established  several 
facts  of  basic  interest  in  the  theoretical  study  of  premixed  flames.  Some  of 
these  have  been  anticipated  before  but  the  present  analysis  appears  to  provide 
the  first  reasonably  sound  mathematical  evidence.  These  results  are  now 
summarized . 

When  the  Lewis  Number  is  equal  to  one,  a flame  that  comes  close  to  a surface 
through  which  there  are  heat  losses  behaves  in  one  of  two  ways.  If  the  heat 
losses  are  small  the  dominant  effect  is  a geometrical  one  which  causes  the 
flame  speed  to  increase.  No  change  of  flame  temperature  is  associated  with  this 
change  of  flame  speed  in  the  sense  that  when  the  heat  losses  are  zero  the  flame 
temperature  is  everywhere  equal  to  the  adiabatic  flame  temperature.  On  the 
other  hand  if  the  heat  losses  are  sufficiently  large  the  flaune  speed  decreases 
as  the  wall  is  approached,  attaining  the  value  zero  at  some  finite  distance 
from  the  wall,  and  this  is  identified  with  quenching.  This  decrease  in  flame 
speed  is  associated  with  a decrease  in  flame  temperature  below  the  adiabatic 
value. 

A flame  of  arbitrary  Lewis  Number  (provided  |l  - l|  = 0(^))  that 
approaches  a shear  flow  can  also  be  quenched,  and  whether  this  occurs  or  not 
depends  on  the  parameter  values.  For  very  large  shear  gradients  quenching  does 
not  occur  for  a finite  band  of  non-positive  Lewis  Numbers,  but  outside  of  this 
band  quenching  does  occur.  In  a neighborhood  of  the  quenching  point  the  flame 
temperature  is  below  the  adiedsatic  value. 

Although  all  of  these  results  are  valid  only  for  Lewis  Numbers  within  O(^) 

U 

of  1,  it  can  not  necessarily  be  argued  that  this  restriction  inevitably 
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•-xcludes  large  numbers  of  combustible  mixtures.  In  practice  activation  energies 
are  only  moderately  large,  with  values  of  9 between  10  and  20  being  fairly 
typical,  and  for  such  modest  numbers  0(1)  values  of  A cover  the  whole 
range  of  Lewis  Numbers  possible  for  a gas  mixture.  Thus  it  is  conceivable  that 
for  steady  premixed  laminar  flames  the  only  results  of  activation  energy 
asymptotics  that  are  of  practical  interest  are  those  for  which  |l  - l|  is 
O(^) . The  stability  results  of  Buckmaster  (1977)  and  Sivashins)ty  (1977) 
support  such  a view. 

In  spite  of  these  remarlts  it  would  be  pranature  at  this  stage  to  discard 
theoretical  results  such  as  those  for  slowly  varying  flames  valid  when  |l  - l| 
is  0(1),  and  one  of  the  results  of  the  present  paper  is  a remar)cably  simple 
physical  interpretation  of  the  general  equation  governing  such  flames.  Thus 
it  is  shown  that  a simple  function  of  the  flame  speed  is  proportional  to  the 
logarithmic  time  derivative  of  a volume  element  of  the  flame,  where  the  sign  of 
the  proportionality  constant  depends  on  the  sign  of  (L-1) . The  role  played 
by  changes  in  the  flame  thic)^ness  as  well  as  the  flame  stretch  destroys  the 
notion  that  useful  insight  into  the  quenching  of  flames  can  be  obtained  by 
considerations  of  stretch  alone;  and  this  objection,  deriving  as  it  does  from 
a rational  analysis  of  the  combustion  equations,  is  quite  independent  of  any 
future  conclusions  that  may  be  reached  about  the  practical  utility  of  slowly 
varying  flame  results. 

There  is  an  important  lesson  to  be  learnt  from  the  history  of  flame  stretch 
of  which  the  present  discussion  forms  part.  Mathematical  studies  are  often 
criticized  for  lacking  physical  content,  and  despite  the  remark  made  in  the 
Introduction  that  equations  such  as  (1.1)  have  a scientific  life  of  their  own 
independent  of  known  physical  facts,  certainly  it  is  wise  to  bear  in  mind  their 
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physical  origin.  But  while  admitting  that  mathematics  without  physics  can  be 
sterile,  it  also  should  be  noted  that  physics  without  mathematics  can  simply 
be  wrong.  The  complexities  of  problems  such  as  those  examined  here  are  too 
great  to  be  unravelled  by  intuitive  reasoning  alone,  so  that  valuable  though 
the  latter  may  be  (the  search  for  eejuation  (A. 19)  was  motivated  by  claims  for 
the  significance  of  stretch) , without  very  strong  experimental  evidence  or 
mathematical  evidence  the  conclusions  of  such  reasoning  should  not  be  ta)(en 
too  seriously. 

We  conclude  with  some  observations  on  possible  extensions  of  this  wor)c. 
The  solution  obtained  in  §2  when  a is  zero  and  there  are  no  heat  losses 
through  the  wall  corresponds  to  a two-dimensional  burner  tip  since  the  wall 
condition  is  then  simply  one  of  symmetry.  It  would  be  of  great  interest  to 
study  such  burner  tips  (and  their  axisymmetric  equivalents)  for  different 
values  of  X to  see  whether  or  not  open  tips  can  be  obtained  under  the 
appropriate  circumstances,  and  such  a study  is  currently  underway.  Also,  as 
noted  in  the  Introduction,  the  present  techniques  can  be  applied  even  when  the 
fluid  mechanics  is  incorporated  in  the  proper  fashion,  and  this  facet  will  be 
explored  in  the  future. 
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dix.  The  Fundamental  Equation  Governing  Slowly  Varying  Flcunes. 


In  recently  publisheKl  work  by  the  present  author  (Buckmaster  1977)  an 
equation  is  deduced  that  governs  the  general  motion  of  a slowly  varying  flame. 
The  equations  that  form  the  basis  of  this  investigation  are  similar  to  those 
of  the  present  paper  but  are  more  general  in  that  the  fluid  mechanics  is 
properly  incorporated  into  the  mathematical  model.  In  discussing  this  work  we 
shall  retain  the  notation  of  the  original  paper  rather  than  that  of  the  present 
^>x>rk  in  order  to  minimize  the  chances  of  confusion.  The  chief  differences  are 
that  x^,  X,  y,  z and  t are  dimensional  quantities,  replaces  T as  the 

nond imensional  temperature,  and  the  flame  is  propagating  from  left  to  right. 

A slowly  varying  flame  is  characterized  by  the  fact  that  its  deformation 


is  described  on  the  length  scale 


and  the  changes  in  this  deformation  on 


the  time  scale  — — . Since  the  analysis  encompasses  large  deformations,  with 
m C 

the  location  of  the^flame  sheet  defined  by 


X = x^(y,z,t)  (A.l) 

the  variable  x^  has  magnitude  • 

P 

The  asymptotic  analysis  must  be  carried  out  in  three  regions:  the  preheat 


zone  of  thickness  ol 1 ; the  flame  sheet  of  thickness  o(— — — -I  ; and 

\ mC  / I 0 roC  ) 

P P 

the  hydrodynamic  regime  where  |x  - x^j  is  hydrodynamic 

P 

scale  the  flame  is  a discontinuity  of  the  kind  described  by  Markstein  (1964) 
across  which  the  temperature  jumps  from  to  • The  solution 

on  this  scale  is  obtained  in  principle  by  solving  Eulers  equations 
of  hydrodynamics  on  each  side  of  the  flame  subject  to  the  requirement  that  the 
mass  and  momentum  fluxes  are  conserved  through  the  discontinuity.  This 
analysis  cannot  be  uncoupled  from  that  of  the  other  regions  and  may  be  thought 
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of  as  specifying  the  gas  velocity  immediately  ahead  of  the  flame  (at 


X » Xj  + 0)  as  a functional  of  x^.  This  velocity  is 

v(y,z,t)  = ■ 

Analysis  of  the  combustion  zone  then  establishes  that  changes  in  are 

governed  by  the  system 


H^lln(H^)  = -B 
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tion  parallel  to  the  x-axis  in  the  sense  that  in  the  preheat  zone  the  temperature 
has  the  form 


~ "P  + exp 


r 4 

±"L 

■ 6^  Xi 


(x-x^) 


Our  discussion  is  concerned  with  che  physical  interpretation  of  these  equations, 
and  their  relation  to  flame  stretch  and  other  )cinematic  concepts. 


The  Kinematics  of  a Flame  Sheet.  We  now  examine  the  flame  on  the  hydrodynamic 
scale  and  calculate  the  nondimens ional  stretch  factor 


' mV  ^ 

m C 

P 

in  terms  of  the  arbitrary  functions  x^,  v.  The  ultimate  goal  is  to  interpret 
the  system  (A. 3)  in  terms  of  kinematic  variables  such  as  S,  so  that  the 
discussion  is  carried  out  in  a Cartesian  frame. 

The  most  transparent  way  to  calculate  S,  though  certainly  not  the  most 
elegant,  is  to  consider  the  deformation  of  an  infinitesimal  parallelogrcun  as 
it  moves  over  the  flame  surface.  To  this  end,  consider  a point  A attached 
to  the  surface  with  coordinates  (x,y,z)  at  time  t so  that 

X = x^(y,z,t)  . 

3Xf 

Consider  also  a nearby  point  B with  coordinates  (x.  + e -r — , y + c,z)  and 

3x^  ' 

a point  C with  coordinates  (x,  + c -r — , y,z  + e) . These  three  points 

I dz 

define  an  infinitesimal  parallelogram  in  the  flame  sheet  of  area  A where 


A = 1 + 


The  velocity  of  A is  labelled  q(y,z,t)  = the  velocities  of 


-»  3q  3q 

B and  C are  respectively  ^ t ^ and  ^ Jz 


Thus  after  a time  6t 


the  coordinates  of  A are  (x^  + y + qjSt,  z + q^At)  and  this  only 

remains  in  the  flame  surface,  as  required,  if 


3t  ^*2  3y  h 3z  * 


3x^  3qj 

The  new  coordinates  of  B are  (x^  + € + qj^6t  + € 6t, 

aq, 

v + f + q6t  + « -T — 6t,  z ♦ q,6t  + e -r — 6t)  with  a similar  result  for  C 

^ ^2  3y  3 
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and  in  this  way  the  new  area  of  the  parallelogram  can  be  calculated,  and  after 
dA  „ 

equating  to  A + — 6t  yields  the  result 


1 ^ 
A dt 
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or,  eliminating  using  (A. 7) 
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A dt 


A = J_  4.  -La.  _L 
Dt  ■ 3t  ^*2  3y  ^^3  3z  ' 

The  velocity  of  the  point  A tangential  to  the  sheet  is  equal  to  the 
tangential  component  of  the  gas  velocity  v so  that 

q - (q.n)n  = v - (vn)n 

where 

9x,  9x, 


(A. 8) 


1, 


n = 
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9z  / J 

is  the  unit  normal  to  the  flame  surface.  Noting  that 

9x 

(q-n)  * 


Li 
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equation  (A. 9)  can  be  used  to  express  q in  terms  of  v.  Thus 
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where  the  operator 


Dt 


can  now  be  written  in  the  form 


_D _3_  ^ _3 L J 1 3 

Dt  3t  '^2  3y  '^3  3z  6^  8y  3y  6^  3z  3z 


(A. 13) 


Note  that  if  and  both  vanish  identically  and  depend  only 

on  y,  the  stretch  factor  may  be  written  as 


A dt  dy 


1 + (x^' 


(l+(xp^)^ 


(A. 14) 


which  differs  significantly  from  the  Karlovitz  Number  (equation  3.3). 

Consider  now  the  variations  in  the  thic)tness  6(y,z,t)  of  the  element  of 
the  flame  as  it  moves  over  the  surface.  In  view  of  (A. 10), 


6 = (A. 15) 
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3x  1^1 

1 4. 
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3y  I 

3z  / _ 

and  moreover 


d6  _ d6 

dt  ~ Dt  ■ 


The  vol’'.roe  V of  the  flame  element,  defined  by 

V = 6A 


changes  according  to 


V dt  A dt  6 dt 


so  thiit  from  (A.  12)  and  (A.  10)  it  follows  that 


i ^ 

V dt 


3^x, 


3y 


(A. 16) 


(A. 17) 


Finally  we  note  that  the  flame  speed  is  defined  as  the  normal  component  of  the 


flame  element  velocity  relative  to  the  gas,  i.e. 


M 

i 
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(A. 18) 


3x^  8x^ 

-r — “V  + V + V — 

9^  ''i  2 9v  3 

Flame  Speed  = (q-n)  - (v-n)  = ^ — (A.l 

* (^f  ^ (^)T 

and  the  thicltness  6 is  inversely  proportional  to  this  quantity. 

Interpretation  of  the  Activation  Energy  Results.  Comparison  of  the  asymptotic 
results  (A. 3)  and  the  )cinematic  expression  (A. 17)  shows  that  the  fundamental 
equation  governing  the  propagation  of  a slowly  varying  flame  can  be  written  in 
the  form 


H^^nCH^)  . 0 


(A. 19) 


where  H may  be  identified  with  the  nondimens ional  flame  speed.  Note  that  the 

effect  of  a change  in  volume  depends  on  the  value  of  the  Lewis  Number  L. 

If  L > 1 then  6 is  negative  so  that  dilatation  > oj  slows  the  flame 

below  the  adiabatic  value  (H^  < 1)  whereas  compression  < oj  speeds 

2 

it  up  (H  >1).  If  the  Lewis  Number  is  less  than  one  this  effect  is  reversed. 
Since  6h  = p /m,  (A. 19)  can  be  written  in  the  form 


thus  establishing  a functional  relationship  between  H and  S. 


(A. 20) 
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